NNLO Corrections to the Higgs Production Cross Section. 
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We discuss the next-to-next-to-leading order (NNLO) corrections to the total cross section for (pseudo-) scalar 
Higgs boson production. The computation is carried out in the effective Lagrangian approach which emerges from 
the standard model by taking the limit mt — > oo where mt denotes the mass of the top quark. 



The Higgs boson (H) has not yet been dis- 
covered. Recently a new central value for the 
top quark mass of m t = 178.0 GeV/c 2 was re- 
ported [1] so that the lower limit on the mass 
of the Higgs from the LEP experiments is now 
to - 117 GeV/c 2 [2]. The search for the Higgs 
boson is a high priority item at the Tevatron and 
at the future LHC. 

In this contribution we discuss H produc- 
tion via the gluon-gluon fusion mechanism. In 
the standard model the H-g-g coupling contains 
quark loops. Since the H — q — q coupling is pro- 
portional to the mass of the quark the top-quark 
loop gives the largest cross section. The calcula- 
tions simplify if one takes the limit m t — > oo and 
uses the heavy quark effective theory (HQET) 
containing a direct Hgg coupling. Studies show 
that this is also an excellent approximation. The 
Lagrangian density is 



with 



C* ff =G H <P li (x)0(x), 



with 



0(x) = ~-G%(x)G a ^(x), 
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whereas pseudo-scalar Higgs (A) production is 
obtained from 
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where $ H (x) and $ A (x) represent the scalar and 
pseudo-scalar fields respectively and rif denotes 
the number of light flavours. Furthermore the 
gluon field strength is given by G% v and the 
quark field is denoted by The factors mul- 
tiplying the operators are chosen in such a way 
that the vertices are normalised to the effective 
coupling constants Gh, Ga and Ga- The latter 
are determined by the top-quark triangular loop 
graph, including all QCD corrections, taken in the 
limit m t — > oo which describes the decay process 
B — * g + g with B = H, A namely 

G B = -2 5 / i a s (tf)G 1 J 2 T B F B (T B ) 
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where a s (/i 2 ) = a s (/i 2 )/47r and fj, r denotes the 
renormalization scale. Further Gf represents the 
Fermi constant and the functions F B are given in 
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Figure 1. The reaction p 1 (Pi) +p 2 {Pi) -» H+'X'. 



the asymptotic limit by 
lim F h (t) = ^- , lim F A (r) = - cot/3. 



(6) 



cot (3 denotes the mixing angle in the Two-Higgs- 
Doublet Model [3]. The coefficient functions Cb 
originate from the corrections to the top-quark 
triangular graph provided one takes the limit 
m t — > oo. They were computed up to order a 2 in 
[4], [5] for the H and in [6] for the A. 

One tree diagram contributing to the reaction 
Pi(Pi) +p 2 (P2) -» H +' X' is shown in Fig. 1. 
On the Born level we have the partonic reaction 



9{Pi) +9{Pi) -> B(p 5 ) • 



(7) 



In NLO we have in addition to the one-loop vir- 
tual corrections to the above reaction the follow- 
ing two-to-two body processes 



9 + 9 
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g + q(q) -» B + q(q) , 



q + q —> B + g . 



(8) 



Such diagrams contain the same basic building 
blocks as in the Drell-Yan production of dileptons 
(see the contribution of van Neerven to these pro- 
ceedings) . The LO and NLO contributions to the 
total cross section for H production were com- 
puted in [7]. The two-to-three body processes 
were computed for H and A production in [8] and 



[9] respectively using hclicity methods. The one- 
loop corrections to the two-to-two body reactions 
above were computed for the H in [10] and the 
A in [11]. These matrix elements were used to 
compute the transverse momentum and rapidity 
distributions of the H up to NLO in [12], [13] and 
the A in [11]. 

In NNLO we receive contributions from the 
two-loop virtual corrections to the Born process 
in Eq. (7) and the one-loop corrections to the 
reactions in Eq. (8). To these contributions 
one has to add the results obtained from sev- 
eral two-to-three body reactions containing glu- 
ons and quarks, e.g. 



9(Pi) + 9{V2) -> g(P3) + g(Pi) + B(p 5 ) 



(9) 



The effective Lagrangian method was also ap- 
plied to obtain the NNLO total cross section for 
H production by the calculation of the two-loop 
corrections to the H-g-g vertex in [14], the soft- 
plus- virtual gluon corrections in [15], [16] and the 
computation of the two-to-three body processes 
in [17], [18], [19], [20]. These calculations were re- 
peated for pseudo-scalar Higgs production in [21], 
[22] and [19], [20]. In the case of A production one 
also has to add the contributions due to interfer- 
ence terms coming from the operators 0\ and O2 
in Eq. (4), (see Fig. lb and Fig. 2b of [21]). 

Let us start with the NNLO virtual corrections. 
The most complicated diagram is the two-loop 
non-planar graph for H(p2—pi) — > g{pi)+g(— P2) 
shown in Fig. 2. The momenta are incoming and 
k\ , &2 label the two loop momenta. 

The basic integral for this diagram has been 
known for a long time. It is 
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where s = (pi — P2) 2 and S n is the angular factor 
in n = 4 + e dimensions. In Drell-Yan production 
of lepton pairs the highest power of momenta in 
the numerator was three which required the result 



{kl.p 2 f Vi23567 = ss e S 2 n 



69 



3 



H(p 2 - pi) 



g{pi -ki- k 2 ) 



9 (Pi - ki 



g(pi) 




g(P2 -h- k 2 ) 



g(p2 - k 2 ) 



prescription for the 7 5 -matrix and the Levi-Civita 
tensor in [25], [26], [27]. 

The evaluation of the two-to-three body phase 
space integrals is possible in suitable frames. 
Since we integrate over the total phase space 
the integrals are Lorentz invariant and therefore 
frame independent. The matrix elements of the 
partonic reactions can be partial fractioned in 
such a way that maximally two factors P^j = 
(pi + Pj) 2 depend on the polar angle 9 and the 
azimuthal angle <j>. Furthermore only one factor 
depends on 9 whereas the other one depends both 
on 9 and (f>. Therefore the following combinations 
show up in the matrix elements 
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Figure 2. The two-loop non-planar crossed graph. 
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However in the Higgs case we need the integrals 
with up to four powers in the numerator such as 
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The virtual corrections were calculated in [14] 
with the program Mincer [23] and we have 
checked them by using the methods in [24]. In 
the pseudoscalar case everyone used the HVBM 



Pi =P 2 j =P 2 m=Pl = , Pl=™ 2 . (13) 

For the first combination it is easy to perform the 
angular integrations exactly in the CM frame of 
the incoming partons since all momenta represent 
massless particles. This requires the result in [24]. 
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where the function dj(9,(j),x) = (1 — cos#)~ fe 
(1 — cos 9 cos x ~ sin#cos</>sinx)~'. Here cos% 
is related to kincmatical variables such as x (= 
m 2 /s, y / (s)=CM energy) and another two inte- 
gration variables z and y, and Fi i2 (a, b, c; t) is the 
hypcrgeometric function. 

The angular integral of the second combination 
is more difficult to compute because one particle 
is massive and the result is an one dimensional in- 
tegral over a hypergeometric function which how- 
ever can be expanded around e. Examples of 
these types of integrals can be found in Appendix 
C of [28]. The last combination is very difficult 
to compute in n dimensions because both factors 
contain the massive particle indicated by p$ . This 
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Figure 3. The total cross section cr tot with all 
channels included plotted as a function of the 
Higgs mass at >/S =14 TeV with fi = m 



Figure 4. The quantity N((j,/no) plotted in the 
range 0.1 < n/no < 10 with ^ = m and m = 100 
GeV/c 2 . 



combination can be avoided if one chooses one of 
the following three frames, namely (1) the CM 
frame of the incoming partons, (2) the CM frame 
of the two outgoing partons indicated by the mo- 
menta P3 and P4, and (3) the CM frame of one of 
the outgoing partons and the H indicated by the 
momenta p4 and p$ respectively 

In (3) we performed one angular integration 
(say (f>) exactly and performed the remaining 
9 integration after expanding the integrands in 
powers of e = n — 4. In all these frames, 
using various Kummer's relations, the hyper- 
gcomctric functions are simplified to the form 
F 1>2 (±e/2,±e/2,l±e/2;f(x,y,z)) which is the 
most suitable for integrations over z and y. The 
higher powers of 1/(1 — zj a+l3e or l/z a+/3e , where 
a > 1, were then reduced by successive integra- 
tion by parts with exact hypergeomctric func- 
tions until we arrived at terms in 1/(1 — z) 1+/3e 
or l/z 1+l3e multiplied by regular functions. This 
required the following identity 

^F li2 ( e /2, £ /2,l + £ /2;/(z)) = 
-Fi, 2 (e/2,e/2,l + e/2; /(*))). 



(15) 



The remaining integrals like Jq dzz 1 l3e f(z) 
and/or ^ dz{\ - z)- 1 -^ f{z) were simplied as 
follows: 

/ dzz- x -^f{z) = 
Jo 

[ dz z' 1 - (3e f{z)+ [ dz z- x - pe f{z) ,(16) 

where 5 « 1. The first term can be evaluated to 
be /(0) ( [-ps] - 1 + In S + [-0e/2] In 2 6 + ■ ■ •) . After 
expanding z~P £ in powers of e in the second term 
the z integration can be performed exactly order- 
by-order in s with non-zero 5. At the end the 8 
dependence cancels in each order in e. Since the 
z integration over the hypergeometric functions 
is nontrivial due to their complicated arguments, 
we expanded them in powers of e prior to the z 
integration, using 

Fi, 2 (e/2,e/2,l+|;/(«)) = l 

+ jU 2 (f(z)) + y (S li2 (/(z)) - Li 3 (/(z))) 

+ ^(Si, 3 (.fW)-S 2 , 2 (/(,)) + Li 4 (/(z))) 

(17) 
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where S„. p (z) is the Nielsen function, with n,p > 
1 and Li„(z) = S„_i i i(^)(see [29]). 

We have repeated the same procedure to per- 
form the remaining y integration. The integrals 
were programmed using FORM [30] by extend- 
ing a program originally used to compute the 
NNLO coefficient functions of the Drell-Yan pro- 
cess in [31], [32]. After performing all these inte- 
grals, we removed all the ultraviolet divergences 
by coupling constant and operator renormalisa- 
tion. The remaining collincar divergences are 
removed by mass factorisation. Then we are 
left with finite partonic cross sections which are 
folded with parton distribution functions to the 
compute hadronic cross section for inclusive H 
and A production. 

Our method differs from the one presented in 
[17], [21] and the approach followed in [18], [22]. 
The authors in the latter references compute the 
total cross section using the Cutkosky [33] tech- 
nique where one- and two-loop Feynman integrals 
are cut in certain ways. These Feynman integrals 
can be computed using various techniques (for 
more details see the references in [18]). The au- 
thors in [17] chose the more conventional method 
which was already used to compute the coefficient 
functions for the Drell-Yan process. However in- 
stead of an exact computation of the 2 — > 2 and 
2^3 body phase space integrals they expanded 
them around x — m 2 / s = 1. Since the coefficient 
functions were known from the Drell-Yan reaction 
to be expressible in a finite number of polyloga- 
rithms like Li„(x), S„ iP (x) and logarithms of the 
types In' x h\ 3 (1 — x), which are all multiplied by 
polynomials in x, one could expand these func- 
tions in the limit x — > 1 and match them with 
the expressions coming from the phase space inte- 
grals. In this way the coefficients in the expansion 
were determined. 

In Fig. 3 we show the total cross section for 
H production in p-p collisions at the LHC. For 
the computation of the effective coupling con- 
stant Gb in Eq. (3) we choose the top quark 
mass m t — 173.4 GeV/c 2 and the Fermi constant 
G F = 1.16639 x 10" 5 GcV~ 2 = 4541.68 pb. We 
used the parton densities in [34]. The cross sec- 
tions for the A are about 9/4 cot 2 (3 times those 
for the H. Next we study their variation with rc- 
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Figure 5. The AT-factors in NLO and NNLO at 
y/S = 14TeV as a function of the Higgs mass us- 
ing the MRST parton density sets; K NNLO solid 
line, K NLO dot-dashed line. 



spect to the scale fi by computing the ratio 



N 
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(18) 



where fiQ = m and /x is varied in the range 
0.1 < ,u/,uo < 10. In Fig. 4 a plot of N is shown 
for m = 100 GeV/c 2 . Here one observes a clear 
improvement while going from LO to NNLO. In 
particular the curve for NNLO is flatter than that 
for NLO. Another way to estimate the reliability 
of this result is to study the rate of convergence 
of the perturbation series, which is represented by 
the A-factor. We choose the following definitions 
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In Fig. 5 one observes that both if -factors vary 
slowly as m increases. Moreover there is consider- 
able improvement in the rate of convergence if one 
goes from NLO to NNLO. At m = 100 GeV/c 2 
we have A NLO ~ 1.74 whereas K NNLO ~ 2.00. 
Still the corrections for Higgs boson production 
are larger than those obtained from Z-boson pro- 
duction at the LHC where one gets A NLO - 1.22 
and K NNLO ~ 1.16 (see [32]). 
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